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ABSTRACT 

We use deep Hubble Space Telescope (HST) imaging of the Frontier Fields (FF) to accurately measure 
the galaxy rest-frame ultraviolet luminosity function (UV LF) in the redshift range z ^ 6 — 8. We 
combine observations in three lensing clusters A2744, MACS0416, MACS0717 and their associated 
parallels fields to select high-redshift dropout candidates. We use the latest lensing models to estimate 
the flux magnification and the effective survey volume in combination with completeness simulations 
performed in the source plane. We report the detection of 227 galaxy candidates at z = 6 — 7 and 25 
candidates at z ^ 8. While the total survey area is about 4 arcmin^ in each parallel field, it drops to 
about 0.6 to 1 arcmin^ in the cluster core fields because of the strong lensing. We compute the UV 
luminosity function at z ^ 7 using the combined galaxy sample and perform Monte Carlo simulations 
to determine the best fit Schechter parameters. We are able to reliably constrain the LF down to 
an absolute magnitude of Muy = —15.25, which corresponds to 0.005L*. More importantly, we find 
that the faint-end slope remains steep down to this magnitude limit with a = —2.041q ]^ 7 . We find 
a characteristic magnitude of M* = — 2Q.89tn!72 Lo g(0^) =—3.54tQ;45. O ur results confirm the 

most recent results in deep blank fields (jFinkelstein et al.ll20l4 iBouwens et al.ll^l5b[) but extend the 
LF measurements more than two magnitudes deeper. The UV LF at z ~ 8 is not very well constrained 
below Muy = —18 due to the small number statistics and incompleteness uncertainties. To assess 
the contribution of galaxies to cosmic reionization we derive the UV luminosity density at z ^ 7 by 
integrating the UV LF down to an observational limit of Muv = —15. We show that our determination 
of Log(pt/v)=26.2 zb 0.13 (erg s“^ Hz“^ Mpc“^) can be sufficient to maintain reionization with an 
escape fraction of ionizing radiation of fesc = 10 — 15%. Future HFF observations will certainly 
improve the constraints on the UV LF at the epoch of reionization, paving the way to more ambitious 
programs using cosmic telescopes with the next generation of large aperture telescopes such as the 
James Webb Space Telescope and the European Extremely Large Telescope. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: luminosity function — 
gravitational lensing: strong 
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1. INTRODUCTION 

One of the most important challenges in observational 
cosmology is the identification of the sources responsible 
for cosmic reionization. Shortly after the Big Bang the 
Universe was completely neutral following the recombi¬ 
nation of hydrogen atoms, until the first sources started 
to reionize the neutral gas in their surroundings. Several 
observational results have now narrowed down the period 
of reionization to the redshift interval 6 < z < 12. Ob¬ 
servations of the Gunn-Peterson effect in the absorption 
spectra of quasars and gamma-ray bursts (GR B) indicate 
that the Universe was most ly ionized by z ^ 6 (|Fan et al.l 
l2006l : [Chornock et al.ll201^ . The sensitivity of Lya emis¬ 
sion to neutral gas is also used to probe the ioniza¬ 
tion state of the intergalactic medium (IGM) at z > 6 . 
In particular, the prevalence of Lya emitters (LAEs) 
among continuum selected galaxies (LBGs) appears to 
drop very rapidly at z > 6.5, which s uggests an iiicrease 
in the fraction o f neutral hydrogen (iStark et al.l 120101: 
Treu et all 120131 ISchenker et al.ll2Q14at IPentericci et al l 


2014f) . TTe decline in the LAEs fraction also favors a 


patchy rather than a smooth reionization process. The 
optical depth of Thomson scattering to the cosmic mi¬ 
crowave background (GMB) recently reported by the 
Planck collaboration ( Planck Collaboration et al.H2015ll 


implies a redshift of instantaneous reionization around 
Zr = 8 . 8 ^i[' 2 , significantly lower than earlier determina¬ 
tions of Zr = 10.6 ± 1.1 by the Wilkinson Mic rowave 
Anisotropy Probe (WMAP, iBennett et al.l 120131) . De¬ 
spite these major advances, large uncertainties remain 
regarding the main sources that drive the reionization 
process. 

Early star-forming galaxies are now thought to be 


power (e.g. IRobertson et al.ll2014 

: iMitcheli-Wvnne et ai.l 

20151: IDuncan & Conselicdl2015l 

Finkelstein et al.l 12014 

Mason et al.l 120151). Deep imaging campaigns of blank 


fields with the Hubble Space Telescope and ground- 
based instrumentation have made important inroads in 
constraining the galaxy ultraviole t luminosity function 


Bunker et al.1 120101: 

Oesch et al.l 20121: McLure et al.l 

20131: Schmidt et al 

12014 IBouwens et al.l I2015bl: 

Finkelstein et al.l 120141). which, in turn, encodes im- 


portant information about the cosmic star formation 
(|Bouwens et al.l l2015bl : iRobertson et al.l I2015D . Indeed, 
the rest-frame UV radiation traces the recent star for¬ 
mation averaged over the past hundreds Myr. Therefore, 
the integration of the UV LF provides the luminosity 
density of galaxies at a given redshift, which can be 
translated to the ionizing radiation, assuming a certain 
star-formation history. The UV luminosity density is 
sensitive to two main parameters: (i) the faint-end slope 
of the LF, and (ii) the integration limit at the faint end. 

The faint-end slope of the LF evolves with redshift, 
and gets as steep as a ~ —2 at z ~ 7. For compari¬ 
son, the faint-end slope at z < 0.5 slope is very shallow 
(a ^ —1.3), which indicates a decreasing contribution 
of faint star-forming galaxies to t he total star forma¬ 
tion density towards low er redshift (|Arnouts et al.l[2005l : 
ISchiminovich et aTll2005l) . The predicted evolution of the 
dark matter halo mass function based on cosmological 
simulations predicts an even steeper slope of a ~ 2.3 of 


the UV LF at high redshift (|Jaacks et al.l[20l^ . The 
deepest HST observations of blank fields, such as the 
XDF, put constraints on the faint-end slope down to an 
absolute magnitude of Mijv ~ —17.5 AB. However, the 
ability of galaxies to reionize the Universe relies on the 
extent of the steep faint-end slope down to lower lumi¬ 
nosities, typically around O.OOIL* at z ^ S While 
cosmological simulations point to a halo mass limit of 
lO^MfT^ for ea.rlv galaxy formation (iJaacks et al.l 120121 : 
iKimm fc Cenl 12014 IWise et al.ll2014D . the depth of cur¬ 
rent observations however limits the exploration of the 
LF to galaxies brighter than ~0.1L*. One particularly ef¬ 
ficient way to push the limits of current facilities is to take 
advantage of the gravitational lensing offered by massive 
galaxy clusters, which act as natural telescopes by boost¬ 
ing t he flux of background sources (iKneib fc Nataraiaiil 
[ 2 ^ . 

Since the dis covery of the first g iant arcs created by 
strong lensing (iSoucail et al.l I19CT . cosmic lenses have 
been successfully used to detect intrinsically faint back¬ 
ground sources and perform spatially detailed analysis of 
distant galaxies. Recently, using HST imaging of 2 5 X- 
ray-sele cted clusters in the Abell jAbdle^^lJ 19891) and 
MAC S (lEbeling et al.ll20niL[2?)(T7ll20inHMann fc Ebeline 


2012 ) catalogs, the CLASH program (iPostman et al. 


20121 ) has made important progress in the character¬ 
ization of the lensing properties of clusters. Multi¬ 
wavelength observations has enabled the measurement 
of the total cluster mass with a precision of 10 % and the 


detec tion of some of th e most distant galaxi e s at z > 7 
(e.g. f^rin et al.l 120111: iBradlev et al.l 1201^ iCoe et aH 
120131) . However, two limitations prevented such pro¬ 
gram from exploring the faint-end of the LF at the epoch 
of reionization. First, the relatively shallow data com¬ 
pared to HUDF, restricted the survey to the brighter 
end part of the LF, even in the case of high magnifica- 
tions. In this sens e, wide area surveys, such as BoRG 
(|Trenti et al.ll20lll) . are indeed well suited to help con¬ 
strain the brighter p art of the LF (IBradlev et al.l 1^0121 : 

I Schmidt et aljl2M^ . Second, the lack of very accurate 
lensing models for some of these clusters, mostly due to 
the low number of multiple images available to constrain 
the mass distribution, also thwarts the construction of a 
reliable LF. 

The Hubble Frontier Fields (HFF) project aims at 
overcoming these two major limitations by obtaining 
deep multi-wavelength observations of six massive galaxy 
clusters that act as cosmic lenses. The HFF include 
the deepest optical and near-infrared observations of 
lensing clusters using HST director discretionary time, 
which are complemented by a wealth of data including 
ALMA, Spitzer, Chandra, XMM, VLA, as well as 
HST and gro und-based imaging and spectroscopic 


follow-up fe.ff. lOwers et al.l 12011 

; lEbeling et al.l 12014 

Rawle et al.l 2015: Karman et al. 

120151: Richard et al.l 

20151: iMedezinski et al.l 2015: 

iRodnev et al.l 2015: 

Zitrin et al.l 20151: Schirmer et al.l 12015: Grille et al.l 

20151: iTreu et al.l 120151: lOgrean et al.l I2015D. Based 


on the full HFF dataset of the first cluster A2744, 
lAtek et ahl (I2015D presented the first cons traints on the 
UV LF at z ^ 7 and z ~ 8 (see also llshigaki et al.l 


The characteristic magnitude M^y=-21 AB (e.g. 
IBouwens et al.ll2015bl : lAtek et al.ll2014n) 


































































































































Are Ultra-faint Galaxies at z = 6 — 8 Responsible for Cosmic Reionization ? 


3 


l2015all . The key result was the steep faint-end slope of 
a ^ —2.01 that extends down to Mjjv = —15.5 AB. 
However, large uncertainties arising from small sample 
size, lensing models, and cosmic variance still prevent 
strong conclusions on the total UV luminosity density 
of galaxies at the epoch of reionization. 

In this paper, we combine the complete dataset 
of the three lensing clusters A2744, MACS0416, and 
MACS0717, and their respective parallel fields to search 
for high-redshift dropout galaxies and put stronger con¬ 
straints on the UV LF at z > 6. The paper is organized 
as follows. In Section [21 we describe the observations. 
The sample selection method is described in Section O 
In Section 131 we briefly describe the lensing models and 
the multiple-image identification. The procedure and the 
results of the computation of the UV luminosity function 
are presented in Section[5l A summary is given in Section 
ini Throughout the paper, we adopt a standard ACDM 
cosmology with Hq = 71 km s“^ Mpc“^, Ha = 0.73, and 
Hm = 0.27 to be consistent with previous studies. All 
magnitudes are expressed in the AB system. 


2. HFF OBSERVATIONS 

The HFF clusters and parallel fields were observed by 
HST with three ACS (Advanced Camera for Survey) op¬ 
tical filters (F435W, F606W, F814W) and four WFC3 
(Wide Field Camera 3) near-IR filters (F105W, F125W, 
F140W, F160W). Observations were scheduled in two- 
epoch sequences, obtaining ACS observations of the main 
cluster and WFC3 observations of the parallel field in one 
epoch and swapping instruments in the second epoch. 
We use the high-level science products delivered by the 
Space Telescope Science Institute (STScI) through the 
Mikulski Archive for Space Telescope^ (MAST), which 
include drizzled science and weight images. 

Basic calibrations were performed with the standard 
IRAF procedures CALACS and CALWF3 for ACS and 
WFC3 data, respectively. Here we chose a pixel scale 
of 60 mas pix“^ for both optical and IR drizzled im¬ 
ages. For the ACS bands we used the ’’self calibrated” 
mosaics, which contain additional corrections applied by 
the STScI team to better account for charge transfer in¬ 
efficiency (CTI) effects. Similarly, the WFC3 bands were 
also corrected for a time-variable IR background. We re¬ 
fer the reader to a detailed explanation of the data reduc¬ 
tion performed by STScf^ (Koekemoer et al. in prep). 
For the fields that have been observed prior to the HFF 
program, we combine all the available data using the 
weight maps included in the HFF data release. Table ID 
summarizes the exposure times and the depth achieved 
in each filter for each of the fields. The limiting mag¬ 
nitude in each fitter was calculated using 0.4"diameter 
apertures randomly distributed in the image to sample 
the sky variance before fitting the resulting distribution. 
The quoted depth is given at the 3-cr level. 

3. SAMPLE SELECTION 
3.1. Photometric catalogs 

https://archive.stsci.edu/prepds/frentier/ 

http://www.stsci.edu/hst/campaigns/frentier-fields/ 



Fig. 1. — Throughput curves of the HST/ACS and WFC3 filter 
set used for the HFF observations. 

We constructed the photome tric catalogs in each fiel d 
using the SExtractor software ([Bertin &: Arnoutslll996[l . 
We first matched all the images to the same point spread 
function (PSF) using a model based on t he largest PSF 
derived with TinyTim (iKrist et al.l 1201111 . In order to 
increase the sensitivity to faint sources, we created deep 
images by using inverse variance (IVM) weight maps to 
combine all IR frames for z ~ 7 sources and F125W, 
F140W and F160W for z ^ 8 sources, respectively. This 
deep image is used for source detection in the SExtractor 
dual image mode, while individual images are used for 
photometry, weighted by the individual IVM images. 

One important limitation for the detection of the 
faintest sources in the cluster fields is the contamination 
from the intra cluster light (ICL) in the central region 
of the cluster. This diffuse light, concentrated primar¬ 
ily in the central region of the cluster, is due to the 
brightest cluster galaxy (BCG) or other bright cluster 
members and the tidal stripping of stars from interact¬ 
ing galaxies during the mer ging history of the cluster 
fe.g. iMontes fc Truiillol . In order to mitigate the 
contamination of sources close to the cluster center we 
subtract a median filtered image from the detection one 
with a filter size of ^ 2"x2", while the photometry is per¬ 
formed on the original image. The adopted filter size is a 
tradeoff between the removal of extended bright emission 
and the appearance of artifacts close the bright galaxies 
of the cluster core. We note that the median filtered de¬ 
tection image allows us to detect five more sources on 
average in cluster fields compared to the original images. 
This is also important for the visual inspection of galaxy 
candidates as the background flux is much lower in these 
corrected images. 

In order to estimate the total flux errors, we ran sim¬ 
ulations of galaxies with different sizes and profiles (cf. 
the completeness simulations in Section [0|) and com¬ 
pared the recovered fluxes to the input values. We 
find that the median-filtering approach achieves an un¬ 
certainty of 0.5 mag for the faintest sources (i7i4o ~ 
28 — 29 mag), by using a local background estimate with 
back_size = 6 in SExtractor for photometry. In ad¬ 
dition, the following extraction parameters were set to 
improve the detection of the faintest sources in the field: 
detect_minarea=2 and detect_thresh=l. 5 for the de¬ 
tection and deblend_nthresh=16 for source deblending. 
Two types of photometry are used throughout this work. 
The isophotal magnitude (ISO) is used to compute the 
colors and ensure that the same aperture is used across 
the filter set. The total flux is measured within the Kron 
radius using the AUTO magnitude, which hereafter is 
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TABLE 1 

HST OBSERVATIONS OF THE HFF FIELDS 




A2744 



MACS0416 



MACS0717 


Filter 

#Orbits 

Depths 

Obs Date 

#Orbits 

Depth 

Obs Date 

#Orbits 

Depth 

Obs Date 

WFC3/F160W 

24 

28.3 

Oct/Nov 2013 

24 

29.1 

Jul/Sep 2014 

26 

28.8 

Feb/Mar 2015 

WFC3/F140W 

10 

29.1 

Oct/Nov 2013 

10 

28.8 

Jul/Sep 2014 

12 

28.5 

Feb/Mar 2015 

WFC3/F125W 

12 

28.6 

Oct/Nov 2013 

12 

28.8 

Jul/Sep 2014 

13 

28.6 

Feb/Mar 2015 

WFC3/F105W 

24 

28.6 

Oct/Nov 2013 

24 

29.2 

Jul/Sep 2014 

27 

28.9 

Feb/Mar 2015 

ACS/F814W 

42 

29.4 

Jun/Jul 2014 

50 

29.2 

Jan/Feb 2014 

46 

29.3 

Sep/Nov 2014 

ACS/F606W 

10 

29.4 

Jun/Jul 2014 

13 

29.1 

Jan/Feb 2014 

11 

28.6 

Sep/Nov 2014 

ACS/F435W 

18 

28.8 

Jun/Jul 2014 

21 

30.1 

Jan/Feb 2014 

19 

29.5 

Sep/Nov 2014 



A2744 



MACS0416 



MACS0717 




Parallel 



Parallel 



Parallel 


Filter 

#Orbits 

Depth 

Obs Date 

#Orbits 

Depth 

Obs Date 

#Orbits 

Depth 

Obs Date 

WFC3/F160W 

24 

28.5 

Jun/Jul 2014 

28 

28.5 

Jan/Feb 2014 

26 

28.7 

Sep/Nov 2014 

WFC3/F140W 

10 

28.2 

Jun/Jul 2014 

12 

28.6 

Jan/Feb 2014 

10 

28.3 

Sep/Nov 2014 

WFC3/F125W 

12 

28.34 

Jun/Jlu 2014 

12 

28.6 

Jan/Feb 2014 

14 

28.5 

Sep/Nov 2014 

WFC3/F105W 

24 

28.6 

Jun/Jul 2014 

28 

28.9 

Jan/Feb 2014 

24 

28.6 

Sep/Nov 2014 

ACS/F814W 

42 

29.0 

Oct/Nov 2013 

42 

28.9 

Jul/Sep 2014 

43 

28.7 

Feb/Mar 2015 

ACS/F606W 

16 

28.9 

Oct/Nov 2013 

10 

28.7 

Jul/Sep 2014 

16 

28.4 

Feb/Mar 2015 

ACS/F435W 

28 

29.7 

Oct/Nov 2013 

18 

29.1 

Jul/Sep 2014 

26 

29.1 

Feb/Mar 2015 


®The depth of the images are S-cr magnitude limits measured in a 0.4”aperture. 


used as the total magnitude. We modified the magnitude 
errors to account for pixel-t o-pixel noise cor r elatio ns in 
the drizzled images following iCasertano et al.l (|200Clll . Fi¬ 
nally, the individual catalogs were matched into a master 
photometric catalog and cleaned from spurious sources. 

3.2. High-redshift Dropout Selection 

We adopted the Lyman break se lection technique 
(|Steidel et al.l 119961 : iGiavalisc^ I2002D to detect high- 
redshift galaxy candidates. The selection is based on 
color-color criteria to sample the UV continuum dropout 
of the IGM absorption blueward of Lya caused by the in¬ 
tervening hydrogen along the line of sight, and minimize 
contamination b y low-redshift red objects at the same 
time. Following (|Atek et al.ll2015ll . we use the following 
criteria to select z = 6 — 7 galaxies: 

(.^814—bios) > 1-0 

(^814 —bios) > 0.6-h 2.0(bioS“<^125) (1) 

(bios~>^i2s) < 0.8 

In addition, we require all sources to be detected in 
the deep IR image and in at least two IR bands with 
4-cr significance or higher. We also reject any galaxy 
that shows up at a significant (at l.Scr) level in the deep 
optical combination of B43S + Veoe images. To preserve 
the Lyman break criterion in the case of non detection 
in the /814 image, we assign a 2 -ct limiting magnitude to 
this band. Therefore, we select only galaxies with Fios 
at least 1 mag brighter than the /814 depth. Similarly, 
we select z ~ 8 galaxies that satisfy: 


TABLE 2 

Number of galaxy candidates in each field 


Field 

z = 6- 7 

z = 8 

A2744 

45 

7 

MACS0416 

33 

3 

MACS0717 

41 

3 

A2744 par 

44 

3 

MACS0416 par 

33 

5 

MACS0717 par 

31 

5 


sources (detected in the IR) accepted in our sample. In 
Figure!^ we show the result of our selection procedure. 
The green circles show the location of all high-redshift 
galaxies identified in this work in the color-color dia¬ 
gram. The dotted and solid lines represent the color evo¬ 
lution of low-redshift elliptical galaxies a nd high-redshift 
starbu rsts, respectively, constru cted from iGoleman et ^ 
(|1980ll and iKinnev et al.l (|1996ll templates. An attenua¬ 
tion of Ay = 1,2,3 is applied to the blue, orange, and 
red curves, respectively. The shaded region represent our 
adopted selection window, which was chosen to minimize 
contamination from low-redshift interlopers and red ob¬ 
jects such as cool stars represented by magenta points. 
Our final sample, combining cluster and parallel fields, 
contains 227 galaxies at z ^ 6 — 7, according to the selec¬ 
tion based on Eq. [U and 25 galaxies at z ^ 8 based on 
the selection of Eq. [2] We find roughly the same num¬ 
ber of candidates in the cluster and the parallel fields. 
While the survey area is smaller in the lensed field, the 
magnification bias balances the number density, allowing 
the detection of much fainter galaxies than in the parallel 
fields. 


(bl05~>^125) > 0.5 

(bl05~>^125) > 0.3 -I- 1.6( J125 —F7i4o) (2) 
(</l25“^^14o) < 0.5 

We require a 4-cr detection in a deep Ji25 + I?i40 + -ffi60 
image and the detection in the stacked deep optical image 
(four ACS filters) to be less than 1.5cr. The deep optical 
image is about one magnitude deeper than the faintest 


3.3. Sources of Contamination 

We now discuss the main sources of contamination for 
our high-redshift sample of galaxies. A first possibil¬ 
ity is spurious sources produced by detector artifacts, 
diffraction features, or photometric errors scattering into 
the color selection space. We have visually inspected all 
of our candidates to identify such contaminants. Most 
of the spurious sources were already cleaned using the 
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Fig. 2. — The Color-color selection windows (represented by the shaded regions) of high-z dropout candidates. The left and right panels 
show the selection of z ~ 7 and z ~ 8 galaxies, respectively. The top panels are for cluster fields and bottom panels for parallel fields. 
The green circles (with associated I-ct uncertainties) represent a compilation of all galaxies satisfying our selection criteria and included in 
our z > 6 samples. The d otted line show the color track, i.e. the redshift-evolution of colors, o f low-redshift ellipti cal galaxies generated 
from IColeman et al.l 11198011 templates, whereas solid lines show starburst galaxies generated from IKinnev et al.l 1119991 templates. The color 
code, from blue to red, illustrates the impact of extinction in steps of Ay = 1. Finally, the magenta points denotes the color track of stars 
generated from IChabrier et al.l 1120091 templates. 


weight maps that exclude the frame edges and the IR 
blobs identified by STScI calibrations. We find that the 
remaining artifacts are mostly diffraction spikes of bright 
stars and few very bright cluster galaxies in the field. 
Fake sources due to photometric noise are also minimal 
since we require the detection in at least two different 
bands. 

The second potential source of contamination are low- 
redshift galaxies that show similar colors to high-redshift 
candidates. According to our stringent spectral break 
criteria, such sources need to show an extremely red 
continuum or a large Blamer/4000A break. Dust ob¬ 
scured galaxies would be excluded by our selection be¬ 
cause they would exhibit relatively red colors redward of 
the break. The remaining galaxies that can enter the se¬ 
lection need to have a large Balmer break and relatively 
blue continuum at longer wavelengths. The existence of 
these peculiar objects have been discussed in iHaves et al.l 
(j2012t) where an SED of a young burst superimposed on 
an old stellar population with an extreme 4000A break 
can mimic the Lyman break colors. Such objects should 
be very rare and represent only a minor contamination. 
Moreover, in the case of a significant contamination by 


these interlopers, we expect to detect the blue continuum 
by stacking the optical images of our candidates, which 
is not the case. 

Alternatively, the Lyman break colors can also be 
mimicked by extremely strong emission lines in z = 
1 — 3 star-forming galax ies fe.g.. lAtek et al.l[201lL l2014at 
Ivan der Wei et ahl 120111) . In this class of galaxies the 
contribution of the emission line flux to the total broad¬ 
band flux is about 2 5% on average and can reach 
85% (lAtek et al.ll201lll and can have important impli¬ 
cations, not only for the high-redshift galaxy selection, 
but als o on the age and stelly ma s s derived from SED 
fitting (iSchaerer fc de BarrosI 120091: iWilkins et al.l 120131: 


Schenker et al.l l2014at iPenin et alTl2015t iHuang et alJ 

lT our case, the contamination of one band 


can lead to an artificial Lyman break, because the faint 
continuum remai ns undetected blueward of the break. 
lAtek et all (120111) estimate that the optical data should 
be about 1 magnitude deeper than the detection band to 
be able to rule out such interlopers. Our criteria impose 
a minimum break of 0.8 mag between the object flux 
and the limiting magnitude of the optical band. More¬ 
over, such a contamination should be even smaller be- 
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Fig. 3.— The redshift selection function for all the fields at 2 ~ 7 
and 2 : ~ 8. The curves are the result of our completeness simula¬ 
tions in recovering the input sources marginalized over the redshift, 
while the shaded histograms represent the actual selected galaxies 
in our sample. The median redshift of the our 2 ~ 8 is significantly 
lower than expected due to the depth of the observations and the 
redshift evolution of the luminosity function. 


cause our stacked candidate images, which are at least 
1.3 mag deeper than the candidate’s flux, show no sig¬ 
nificant detection in the optical bands shortward of the 
break. 

In general, observationally, estimating the contamina¬ 
tion rate of such sources proves very difficult at z > 6 be¬ 
cause it hinges on large spectroscopic follow-up programs 
that aim at detecting the redshifted Lya emission line in 
these galaxies. In addition to the large amount of tele¬ 
scope time needed to reach the required depth, the ab¬ 
sence of Lya emission does not exclude the high-redshift 
nature of the associated objects. Indeed, the increasing 
neutral hydrogen fra ction at z > 6 easily absorbs and 
diffuses Lyg photonsjiAtek et al.ll2009l: iStark et al.ll201Cll : 
iCariiana et al.l[2C)l^ ISchenker et al.ll2014b[ l. ThCTefore. 
one must rely on simulated colors based on galaxy spec¬ 
tral templates to estimate the contamination rate of 
low redshift galaxies, which thus has been found to be 
less than 10 % in mo st studies (e.g. lOesch et al.l 120101 : 
iBouwens et aLll^OlSbU . 

Another source of contamination to consider are low- 
mass stars. As we can see in Fig. [2] a significant fraction 
of these stars (shown in magenta) can have similar colors 
to high redshift galaxies. However, we have excluded any 
source that has a SExtractor stellarity parameter greater 
than 0.8 to minimize point-like objects in our sample. 
Also, the number density of cool stars ranges from 0.02 
to 0.05 arcmin^ derived from observat ions with similar 
or better depths (jBouwens et ahll^lbH l . This translates 
to about 1 contaminant in the entire z 7 and z ^ 8 
galaxy samples. We also consider transient events as 
possible contaminants. Because observations are taken 
at two different epochs (cf. Table S]), a supernova ex¬ 
plosion may appear in IR images and not in the optical 
ones and could be selected as a dropout candidates. Sim¬ 
ilar colors can also be obtained in the case where the IR 
data were taken hrst and the supernova faded until the 
optical images. Again, different reasons point toward 
a negligible contamination level from transient sources. 


Such objects would show point-like profiles that would be 
excluded by the stellarity criterion explained above and 
our additional visual check for unresolved sources. The 
detection rate of supernovae is also very small and, so 
far, none of SN detecti ons in the dedicated search pro¬ 
gram of the HFF (e.g., iR.odnev et al.l[2?)T5ll were inad¬ 
vert entl y selected as high-z galaxies in the different HFF 
studies (lAtek et al]l20i4bl 120151 : iFinkelstein et al.ll20l4 
llshigaki et al.ll20i5^ 


4. HFF CLUSTER MASS MODELS 


In order to exploit the full potential of the HFF clus¬ 
ter lenses, we first need a robust model describing the 
total mass distribution and the lensing properties. In an 
effort to provide the community with all the required 
lensing maps to interpret background source observa- 


to the HFF fe.e'..IBrar 

ac et al.ll2005l 

IMerten et al.ll20III: 

Johnson et al.l I20I4 

Richard et al.l 

120141: IGrillo et al.l 

2014 ICoe et al .11201 5 

). With the availability of new 


deep HST, Spitzer, and spectroscopic observations, the 
models were signihcantly improved thanks in particular 


(IDieero et al.l 1201^ 

t IGrillo et al.l 120151: IZitrin et al.l 120141: 

Wantr et al. 11201,51: 

Jauzac et al.ll201,5ll20l4llshie-aki et al.l 

20I5all. Most notablv. using a set of ^ 180 and 200 


multiple ima g es in A2744 and MACS04I6, respectively, 
iJauzac et al.l l|20I5l I20I4I1 reconstructed the projected 
cluster mass down to a precision of ^ 1%, which repre¬ 
sents a significant improvement over pre-HFF mass mod¬ 
els. 

In the present study, we use the most recent cluster 
models based on full-depth HFF observations and con¬ 
structed by the CATS (Clusters As Telescopes) team. 
The models developed for the first two HFF clusters, 
namely A2744 and MACS04I6, are already available on 
the MAST archive. However, STScI has started a new 
mass mapping initiative so all the teams provide HFF 
mass models of both clusters taking advantage of the 
full depth of the HFF data, but using the same inputs. 
These ‘unified’ models will be submitted to STScI as part 
of our answer to the HFF call for the community to pro¬ 
vide updated lensing map0. The mass modeling of the 
third cluster MACS0717 will be presented in details in 
Limousin et al. (2015). 

The mass reconstruction of each of the clusters is ex¬ 
plained in the publications listed above. Here we briefly 
describe the key points of the mass modeling procedure. 


The CATS HFF m ass models a r e built usi n g the 
LenstooP^ so ftware (|Knei B IT^ l.billo et al.l l2?M 
I Julio k. Kneibl[200^ . For the strong-lensing analysis, we 
are using a parametric approach that consists of model¬ 
ing the cluster mass distribution using both cluster-scale 
and galaxy-scale halos. Galaxy-scale halos are important 
in the mass modeling, as multiple image configurations 
are impacted by their location. Indeed if a multiple im¬ 
age is close to a cluster galaxy, then the distortion created 
by the lensing effect will be a combination of both the 
cluster potential itself, as well as a galaxy-galaxy lensing 
effect due to the potential of the cluster member.The po¬ 
tentials describing these components are modeled using 


https: //eirchive. stsci. edu/prepds/frontier/lensmodels/ 
http: //projets . leun. fr/projects/lenstool/wiki 
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pseudo-isothermal ellip tical mass distribution (PIEMD, 
lEliasdottir et al.ll2007li . 

For MACSJ0416, we used 149 of the most secure mul¬ 
tiple images (over the 194 identified) to constrain the 
mass model. Our best fit mass model of the cluster inner 
core consists of two cluster-scale halos, well-aligned with 
the light peaks from the two brightest cluster galaxies 
(BCGs), plus 98 galaxy-scale halos, corresponding to the 
brightest cluster galaxies. For Abell 2744, we used 154 
most secure multiple images (over the 181 identified). 
Our best-fit mass model consists of two cluster-scale ha¬ 
los to describe the dark matter distribution on large- 
scale, combined with 733 galaxy-scale halos, describing 
the distribution of cluster galaxies. Both models deliv¬ 
ered really good errors in the predictions of multiple im¬ 
age positions, an rms of 0.68" and 0.79" for MACSJ0416 
and Abell 2744 respectively. For MACS J0717, we used 
the 140 most secure multiple images (over 163 identified) 
to constrain a mass model composed of 4 large scale DM 
haloes, plus 92 galaxy scales haloes, corresponding to the 
brightest cluster galaxies. More details of the SL analysis 
will be given in Limousin et al. (2015) 



0 1 2 3 4 5 

g [mags] 


Fig. 4.— The cumulative surface area in the source plane at z ~ 7 
as a function of the amplification factor (in magnitudes) derived 
from the mass modeling of the three HFF clusters. Uncertainties 
in the surface area are also shown at the 1-cr level . 


4.1. Multiple Images 

In addition to flux magnification, strong lensing pro¬ 
duces multiple images of the same background galaxy. 
Therefore, we need to identify multiple images to be re¬ 
moved from the galaxy number counts before computing 
the UV LF. For each galaxy, we predict the position of 
potential counter images by using the mass model and 
Lenstool to project its position into the source plane be¬ 
fore lensing back the source into the image plane, where 
Lenstool predicts the position of all the multiple images. 
In the vicinity of these positions, we look for dropout 
sources that show similar colors and photometric red- 
shifts. In the case of well resolved sources, we also ver¬ 
ify that they have similar morphological and geometrical 
symmetries constrained by the lensing model. In A2744 
we find three systems at z ^ 7, with a total of nine mult i- 
ple images, which were discussed in lAtek et al.l (|2014bD . 
In MACS0416, we identify seven systems at z ^ 7. 


In the case of MAGS0717, the multiple image region 
extends beyond the WFG3 field of view. Therefore, we 
expect each of the galaxy candidates to have counter¬ 
images with similar flux/magnification ratios given the 
lensing geometry of this cluster. The uncertainty on 
the image position predicted from the lensing model of 
MACS0717 is about rms = 1.9 arcsec, much larger than 
in the two other clusters (cf. Limousin et al. 2015). 
Moreover, the shape of the critical line is not as well 
constrained, which makes it harder to identify multiple- 
image systems. From the Lenstool model, we estimate 
the multiplicity in MAGS0717 to be three on average in 
the WFC3 field of view. Therefore, we divide the num¬ 
ber of galaxy candidate by this multiplicity to obtain the 
final number counts used for the LF calculation. 

5. THE GALAXY UV LUMINOSITY FUNCTION 

We now turn to computing the galaxy UV luminosity 
function at z 7 and z ^ 8 by combining the high- 
z candidate samples from all the fields. For each field 
we calculate the completeness function and the effective 
survey volume to derive the LF. While in the parallel 
fields we proceed with standard completeness simulations 
widely used in analyzing blank field surveys, we need to 
take into account the lensing effects in the cluster fields. 

The primarily goal of the HFF program is to extend 
the current limits of deep galaxy surveys by using strong 
lensing to magnify intrinsically faint sources behind the 
galaxy clusters. In parallel, the survey area in the source 
plane is significantly reduced for high magnifications. 
Consequently, the efficiency of a given lensing cluster in 
probing the high-redshift Universe, which can be quan¬ 
tified by the number of magnified sources discovered, is 
the result of a trade off between the amplification factor 
fi and the source plane area a. We show in Fig. |4] the 
cumulative survey area as a function of the amplification 
factor. 

While the survey area in the image plane corresponds 
to the WFC3 field of view, i.e. about 4.7 arcmin^, we 
can see the total survey are a is reduced to ab out 0.6 
to 1 arcmin^ in cluster fields. iWonsr et~aI1 (|2ni2[) define 
the cluster magnification power as the cross section for 
magnifying a source above min imum threshold of /r = 
3 (see also [Richard et all I2014D . Although the precise 
value of the threshold is somewhat arbitrary, the number 
of high-redshift detections in our galaxy samples peaks 
around ^ = 3, or /i ^ 1.2 mag (cf. the top panel of Fig. 
El For this typical value, the total survey area drops to 
values of around 0.2 to 0.5 arcmin^ for the three clusters 
of this study. 

The parallel fields are located about 6 arcmin away 
from the cluster lenses with a small yet non negligible 
magnification. It is important to assess the lensing ef¬ 
fects on the UV LF, both on the magnitude and on the 
survey volume for these regions as well. Among the dif- 
ferent lensing model s available in the HFF project, only 
iMerten et al.l (120111) provide a wide-held magnihcation 
map for the three clusters that covers their flanked Helds. 
We used the magnihcation maps at z 9 that have a 
resolution of 25 arcsec pix“^ to estimate the hux ampli- 
hcation of galaxies. Since no dehection map is available, 
the volume reduction is estimated by dividing the area 
by the amplihcation map. The harmonic mean of the 
amplihcation factor ranges from 1.11 to 1.23 in the three 
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fields, with typical errors of about 10%. The inclusion 
of the lensing effects introduce only small changes in the 
UV LF, within the error bars, basically shifting Mjjv to 
slightly fainter values and (j) to higher values. 

We now describe how we estimate the effective survey 
volume for each cluster by combining the source plane 
area with the redshift selection function and the recovery 
rate of simulated galaxies as a function of different galaxy 
and lensing parameters. 

5.1. Completeness Simulations 








M.b 


Fig. 5. — The effective survey volume as a function of absolute 
UV magnitude for the redshift 2 ~ 7 sample. Each curve is based 
on our completeness simulations for each cluster or parallel field. 
The total area for the cluster fields is corrected for lensing effects 
whereas the blank field area is based on the full WFC3 field of 
view (cf. text for a detailed explanation). The blue colored re¬ 
gion represents the 95% confidence intervals of the completeness 
estimate. 


Following lAtek et ^ (|2015D we run extensive Monte 
Carlo simulations to assess the completeness level as 
a function of the intrinsic magnitude. A total of 
10,000 galaxies were simulated for each of the six 
fields. We take the galaxy profile into account by cre¬ 
ating two sam ples of exponential disks and de Vau - 
couleur shapes (iFerguson et al.l[2004 iHathi et al.ll2008ll . 
These galaxy profiles are then distorted by gravitational 
lensing according to our mass model. As for galaxy 
sizes, we adopt a log-normal distribution with a mean 


consistency with previous results (IFereuson et alJ 

2004; 

Bouwens et al.l2004 

: IHathi et al.ll2008l: 

Oesch et all 

201C; 

Grazian et al.ll201lL 

2012l:lHuang et al.l 

2013l:IOesch et al. 

2014f). the distribution is based on 

the sizes derived 


from spectroscopically confirmed LEG samples at z ^ 4 


(|Vanzella et al.l 1200911 . while accounting for a redshift 
evolution of the intrinsic physical size of galaxies with 
a factor of fl + zj~^. HST observations of dropout 
galaxies at z > 6 redshifts also reveal small siz es of 
less than 0.3" (|Mosleh et al.ll20l1: lOno et aT]l2013[l and 
smaller for galaxies fainte r than Mjjv — 21 mag. More 
recently iKawamata et al.l (1201511 measured the size of 
lensed dropout galaxies in the HFF cluster A2744 and 
report similar results. Small sizes ar ound 0.1"h a ve als o 
been observed in lensed galaxies by ICoe et aH (1201311 : 
IZitrin et al.l (1201411 . In addition, based on the results of 
( Huang et al.ll201^ . we adopt a size luminosity relation 
r PC with [3 = 0.25, for our simulat ions (see also 
iMosleh et al.ll20i^ IKawamata et al.ll2015ll . 
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Fig. 6. — The distribution of the amplification factor, expressed 
in magnitudes. Top panel shows the distribution for the z ~ 7 
candidates in each of the lensing clusters. For comparison, the 
bottom panel presents the result of the completeness simulations in 
the three cluster cores, which contain about 80,000 objects each. 

After assigning random redshifts in the range [5.5-7.5], 
we simulate galaxy magnitudes in e ach ff^Tband usin g 
star-forming SED templates from (jKinnev et al.l [199^ . 
In this step, we also assign random intrinsic absolute 
magnitude (in the rest-frame UV) in the range Mjjv =[- 
14,-24] mag. We then include 10 simulated galaxies in the 
actual images of each band for a total of thousand images. 
This is where the cluster and parallel fields are treated 
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TABLE 3 

The best fit 2 ~ 7 Schechter parameters in each individual 

FIELD 


Field 

Mt,v 

[AB mag] 

a. 

loglO <t>* 
[Mpc“®] 

A2744 

-20.92 ±0.64 

-2.03 ±0.13 

-3.56 ±0.45 

MACS0416 

-21.04 ± 1.10 

-2.07 ±0.16 

-3.67 ±0.59 

MACS0717 

-21.13 ± 1.52 

-2.02 ± 0.22 

-3.73 ±0.83 

A2744 par 

-20.99 ± 1.05 

-2.03 ± 0.19 

-3.62 ±0.55 

MACS0416 par 

-20.94 ± 1.19 

-2.02 ± 0.20 

-3.59 ±0.70 

MACS0717 par 

-20.89 ± 1.10 

-2.00 ±0.21 

-3.54 ±0.62 


differently in the simulations. In the cluster fields, galax¬ 
ies are not included directly in the images but are simu¬ 
lated in the source plane. They are lensed into the image 
plane using the corresponding mass model. This way we 
ensure that all the lensing effects including magnification, 
shape distortion, position relative to the critical line are 
fully taken into account. Finally, for all the images, we 
run the same procedure used to select high-z galaxies 
and determine the completeness function, which repre¬ 
sents the fraction of the original galaxies recovered in 
our selection as a function of the parameters described 
above. This completeness function is in turn incorpo¬ 
rated in the computation of the effective volume in each 
magnitude bin following the equation: 

^eff = f f f{z,m,n) dn{n,z) dz (3) 

where fimin is the minimum amplification factor Umin 
required to detect a galaxy with a given apparent mag¬ 
nitude m and Vcom is the comoving volume. f{z, m, /r) is 
the completeness function that depends on the redshift z, 
apparent magnitude m, and amplification factor /x, and 
dn{fj,) is the area element in the source plane, which is 
function of magnification and redshift. 

Figure [5] presents the results of our effective volume 
estimates in each field, marginalized over the intrinsic 
absolute UV magnitude. The extent of each filled region 
represents the 68% confidence intervals. We can clearly 
see the importance of gravitational lensing in extending 
the survey depth to fainter galaxies. While the blank 
fields completeness drop abruptly before Muv = —18 
mag, it becomes shallower in cluster fields and extends 
down to Muv = —15 mag, although at a level of 10% or 
less. 


5.2. The UV LF at z = 6 — 7 

Using the derived effective volume as a function of ab¬ 
solute magnitude we compute the UV luminosity func¬ 
tion following the equation: 


where W and Mi are the number of galaxies and the 
absolute magnitude, respectively, in each magnitude bin. 
Following most of the studies in the literature, we choose 
a bin size of 0.5 mag whereas the magnitude varies from 
one field to another. The individual LF determinations 
in each field are presented in Figure [T] The top pan¬ 
els show the results in the cluster fields while the bot¬ 
tom panels are for parallel fields. The LF extends to 


Mjjv = —18.25 in blank fields, whereas it reaches a 
magnitude of Mjjv = —15.25 in cluster fields thanks 
to the lensing magnification. This gain can already be 
seen in the completeness function results (see Fig. 0 
which reaches fainter magnitudes in cluster fields. Figure 
0 also shows the distribution of the amplification factor 
for the dropout samples selected behind the three lensing 
clusters. The flux amplification ranges essentially from 
/i ^ 1.25 to 75, with a median value of 5.2 , 3.9, and 4.9 
in A2744, MACS0416, and MACS0717, respectively. 

In addition to the sample contamination discussed in 
Sect. 13.31 several sources contribute to the uncertainties 
of the LF data points. The mass model errors, affecting 
the magnification and the survey volume (or the source 
plane area U), are propagated into the LF determina¬ 
tion in the case of cluster fields. We also include Poisson 
errors, and the cosmic variance estima te based on the 
recent results of IRobertson et all (1201511 for both cluster 
and parallel fields. Finally, uncertainties from incom¬ 
pleteness simulations, as seen in Fig. [5] are incorporated 
in the final LF. The bright-end is dominated by cosmic 
variance errors and small number counts, while the faint- 
end error bars reflect mostly the large incompleteness un¬ 
certainties at these faint magnitudes and small statistics. 

We determine the shape of the rest-frame UV 
LF by fitting a Schechter function (jSchechted 1197611 
to our data, which has been extensively used to 
describe the gal axy UV LF acros s a wide red- 


20061 IReddv & Steir 

ell 120091: IBouwens et alJ 2006: 

McLure et al. 120091 

Ouchi et al.l 20091: IWilkins et al.l 


)10t lOesch et al.l l2012l : IWillott et al.l 1201311 . The 
Schechter function can be expressed in terms of absolute 
magnitudes as: 


</.(M) = 


In(lO) 

2.5 




We perform a Schechter fit to the LF in each field 
and find a faint-end slope a between -2.0 and -2.07 (cf. 
Eq. 0. We i nclude in the fit the data points from 
iBoiiwens et al.l (|2ni5bl . black squares) to constrain the 
bright-end of the LF. The faintest bin is reached in 
A2744 at Mjjv = —15.25. There is a hint of a shal¬ 
lower slope at those magnitudes that could be the re¬ 
sult of large uncertainties in the completeness estimate, 
which is typically less then 10% in this bin. As a mat¬ 
ter of fact, a similar decline is observed in MACS0416 
around Muv = —16.25, whereas the slope is constantly 
steep at these magnitudes in A2744. It is clear that the 
constraints on the LF in MACS0717 are not as good 
as in the other cluster fields. This is likely the result 
of larger uncertainties on the lensing model due to the 
complex structure of the galaxy cluster. For instance, 
the model prediction for image position has an rms error 
of ^ 1.9 arcsec, while it is is about 0.7 to 0.8 arcsec in 
MACS0416 and A2744. These kind of deviations in the 
LF data points can also be observed in A2744 when using 
an old model based on pre-HFF observations that simi¬ 
larly yielded large uncertainties. Nonetheless, the overall 
shape of the LF remains consistent with the other clus¬ 
ters’ results. 

We now combine all the LF constraints from the lensed 
and parallel fields to compute the most robust UV LF 
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A2744 



A2744 par 



MACS0416 



MACS0416 par 



MACS0717 



MACS0717 par 



Fig. 7. — UV luminosity function at 2 7 computed individually in each field. The blue circles represent our determination with 1-cr 

uncertainties while the blue solid curve is our best Schechter function fit to the LF. We compare our re sults to previou s literat ure results 
in blank fields. The black sq uares and dashed curv e are from a compilation of HST leg acy fields by [Bour ns et al.l ll2Q15blV We also 
show the LF determination of ibchenker et ^ ll2Q13bl . blue squares and dashed curve) and IJVfcLure et al.l II2013I. green squares and dashed 
curve) derived in the UDF12 field. We also include data points on the bright-end of the LF [magenta squares] from a wide area survey by 
iBowler et alJ ^20111 . 



Fig. 8. — The combined constraints from the clusters and parallel 
fields on the UV luminosity function at 2 ~ 7. The color code for 
the data^oints and the best fit Schechter function are the same as 
in Fig. [3 

at z 7 . The result is shown in Fig. [8l Together 
with our data points, we show the most recent results 
from the literature from the blank fields as described 
in the legend. We ran MCMC simulations with 10® re¬ 
alizations to find the best fit to the LF and estimate 
the uncertainties on the Schechter parameters. We find 
that the UV LF at z ^ 7 has a faint-end slope of 
a = —2.04to'i7, a characteristic M* = —20.89lo;®2 
and Log((/)*)=—3.54 ^q'45 . This is an excellent agreement 


with earlier results presented in lAtek et al.l (1201511 . where 
we computed the UV LF in the HFF cluster A2744. In 
Figureini we show the likelihood analysis of the Schechter 
parameters at z 7, marginalized over two param¬ 
eters at a time. Our results are in good agreement 
with the most recent results reported from the blank 
fields (cf. Table El . In partic ular, the faint-end slope 
is close to the determination of iFinkelstein et ^ (|2014[1 
with a = —2.03^20 iBouwens et al.l (|2015bll with 
a = —2.03±0.14. Our results show only a slightly steeper 
faint-end slope tha n the values reported in the Hubble 
extreme deep field (|Schenker et al.ll2013al : iMcLure et al.l 
1201311 with a = l-SO^o^is and a = l.OOlg'isj respec¬ 
tively, and agree within the reported uncertainties. Also 
using the gravitational lensing o f the first HFF clus¬ 
ter A2744. hshigaki et al.l (I20i5all find a slightly shal¬ 
lower slope of a = —1.94 ^q ®q, which still in agreement 
with our value within the errors.. We also, performed a 
schechter fit while excluding data points from MACS0717 
cluster, which yields very similar parameters: a = —2.03, 
a characteristic M* = —20.86 and Log((/)*)=—3.52. Im¬ 
portantly, we note that the uncertainties on the faint-end 
slope decreases to (Tq, 0 .1. 

Several theoretical models and cosmological simula¬ 
tions produce predictions for the UV LF at high red- 
shift. For instance, our LF results at z ^ 7 are 
in good agreement with th e hydrodynamical simula¬ 
tions of IJaacks et al.l (|2012f) with M* = —20.82 and 
Log((/i*)=—3.74, although they predict a steeper faint- 
end slope of a = —2.30 down to similar low er magnitude 
limit o f Mjjv = —15. Theoretical models bv iMason et ahl 
(|2015ll find a closer slope of a = —1.95 ± 0.17 down to a 
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TABLE 4 

Combined constraints on the UV LF at ^ ■ 


Fig. 9.— Likelihood analysis of the Schechter parameters for the LfV LF at z ~ 7: the faint-end slope a, the characteristic and 

Log(0®*“'’). The density plots show the result of MCMC simulations marginalized over two parameters in each of the three panels.The 
orange curves represent the 68% and 95% confidence contours. 

constraints on the faint-end part of the LF at 2 ; ^ 8. At 
brighter magnitudes the LF is better constrained thanks 
to the addition of the parallel fields and appears in agree- 
m ent with pre v ious r esults in the literature. As discussed 
in lAtek et all (120151! , the redshift 8 galaxy selection in 
A2744 clearly shows an ove rdensity (see also lZheng et al.l 
12014IT^higaki et al.ll2015^ . which translates into a large 
excess in the UV LF. Here we decided not to exclude the 
entire A2744 field in the combined LF to ayoid introduc- 
ing a well known bia s due to cosmic yariance (see also 
llshigaki et al.l[2f)T5bl !. 


Muv 

Logfip) 

iperr 

-20.25 

-3.4184 

0.1576 

-19.75 

-3.0263 

0.1658 

-19.25 

-2.9044 

0.1431 

-18.75 

-2.7418 

0.1332 

-18.25 

-2.3896 

0.1401 

-17.75 

-2.1032 

0.1990 

-17.25 

-1.8201 

0.1940 

-16.75 

-1.7548 

0.1893 

-16.25 

-1.6044 

0.2117 

-15.75 

-1.4012 

0.3123 

-15.25 

-1.4012 

0.3122 


TABLE 5 

Comparison of the best fit 2 ~ 7 Schechter parameters 


Reference 


[AB mag] 


logio d* 

[Mpc-3] 


This work 
^telc^t_^ ((20^) “ 
Jshigak^^t_^^ (201^) “ 
^ouwens_e^^^ (201^) 
^inkelsteii^_e^_^ (20M) 
McLure et al. (2013) 


-20.89:!:^-^^ 

—2 

13 

- 3 . 54 ; 

-20.90to73 

—2 m 
^•^-*--0.28 

- 3 . 55 ; 

-20.45l°i 
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O) 
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+0.57 

-0.57 

+0.10 

-0.20 

+0.27 

-0.17 

+0.41 

-0.26 

+0.28 

-0.45 


^ Using the first HFF cluster A2744 

magnitude limit of Muv = —12, corresponding to a halo 
mass of 10® Mq. Similarly, iKimm fc CenI (I2014D find a 
theoretical faint magnitude limit of Mjjv = —13 with 
a faint-end slope of a = —1.9. Semi-analytical models 
of iDayal et al.l (l2014f ) also find a steep faint-end slope of 
a ~ —2.02 at z ^ 7 in good agreement with our obser- 
yations. 

5.3. The UV LF at z ~ 8 

Regarding the redshift z ~ 8 LF, we followed the same 
procedure used for the z ~ 7 LF. However, in lensing 
clusters, the survey volume at z 8 is much more sig¬ 
nificantly reduced than at z ^ 7. In the extreme case of 
MACS0717, the total survey volume is about 460 Mpc^, 
whereas it reaches 3600 Mpc® in A2744. Therefore, we 
expect lower galaxy number counts at z 8 for high 
magnification values. In total, we detect only two galax¬ 
ies with intrinsic magnitudes fainter than Mjjv = —18, 
where the uncertainties on the LF estimate are very large 
(see Fig. [TOt . Unlike the z ~ 7 LF, we do not have strong 


-3 

-4 

-5 

-6 



Bouwens et al. (2014) 
Schenker et al. (2013) 
McLure et al. (2013) 
Schmidt et al. (2014) 


Fig. 10. — The combined constraints from the clusters and 
parallel fields on the UV luminosity function at 2 ~ 8. The 
red circles repre sent our LF de t ermina tion while the red squares 
are taken from (Schen ker et all Il2013al !. the black squares from 
IBouwens et al.l ll2015blL and the blue squares fromJMcLure_et_^' 
II 2 OI 3 I L The green curve is the unbinned UV LF of ISchmidt et at 


II2014I ) . The best Schechter fits from the same literature results are 
also shown with dotted lines (cf. legend in the inset). 


5.4. Implications for cosmic reionization 

In this study, we reliably extend the UV LF to unprece¬ 
dented depth, and put strong constraints on the faint-end 
slope down to Mjjy = —15.25. The most important re¬ 
sult is that the faint-end slope remains very steep down to 
a luminosity of 0.005L*, the characteristic UV luminosity 
at z ^ 7. Also, thanks to the combination of six fields, we 
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Fig. 11.— The evolution of the UV luminosity density of galaxies with redshift. Our determ ination (blue circle) is b ased on the integration 
of the UV LF down to Mjjv = —15. The brown squares repres ent the most recent resu lts of IBouwens et all Il2ni5bl 'l who integrate the UV 
LF down to Muy = —17. The orange triangles are the result of IIFinkelstein et al.l2014ll who integrate to conservat ive l imit of Mjjy = —18 . 
The red diamond and the green asterisk are the UV luminosity densities derived at z ~ 9 by IMcLeod et al.l II2015I1 and lOesch et al.l II2014I1 , 
respectively. Both studies integrate the UV LF down to an observational limit of Muy = —17.7, which correspond to a star formation 
rate of SFR;;^ = 0.7 Mq yr“^. All the error bars indicate 1-cr uncertainties. The shaded grey region show denote the UV luminosity 
density based on the ionizing emissivity IIMadau et al.l[T999ll required to maintain the IGM ionized at a certain redshift for and escape 
fraction of ionizing radiation of /esc ~ 20% and three different values for the clumping factor Chu = 1, 3, and 10. The red shaded region 
shows the 68% confide nce interval for the evolution of the galaxy UV luminosity density based on constraints on the ionizing emissivity of 
IBouwens et al.l (I2015all . The green hatched region corresponds to the UV luminosity density required to maintain the IGM ionized for a 
clumping factor of Ghii = 3 and escape fraction in the range /esc = 10 — 15%, which encompasses our value at z ~ 7. 


have significantly reduced the uncertainties on the faint- 
end slope to about 5%. With these strong constraints 
in hand, we can now integrate the UV LF to derive the 
UV luminosity density at z ~ 7. The main advantage 
of these observations is the ability to set the integration 
lower limit to Muy = —15. Unlike previous results that 
extrapolate the UV LF to lower magnitude, we use an 
observational constraint to estimate the ultraviolet pho¬ 
ton budget from galaxies. For instance, based on the UV 
LF results in blank fields, IBouwens et al.l (l2015bD report 
a total UV luminosity density of Log(pcry) = 25.98±0.06 
erg s“^ Hz“^ Mpc“^ and iFinkelstein et al.l (l2014f) a value 
of Log(pc/y) = 25.77±0.06 erg s~^ Hz“^ Mpc“^ at z ^ 7. 
Both values are computed down to Muy = —17 AB, 
which is the limiting magnitude of their observations. 
Here we compute the UV luminosity density down to 
Muy = —15 and find Log(pc7v)=26.2 ± 0.13 erg s“^ 
Hz“^ Mpc“^. 

In order to determine whether the UV luminosity pro¬ 
duced by galaxies is sufficient to reionize the IGM, one 
needs to estimate additional parameters, for which ob¬ 
servational constraints remain challenging. First, the 
conversion factor from the UV luminosity to ionizing ra¬ 
diation ^ion (ergs“^ Hz) depends on the star-formation 
history of galaxies. The value of ^ion is generally con¬ 
strained using stellar population models of early galax¬ 


ies and the observed UV slope 0 oi z > 6 galax- 


ies (iBolton & Haehnelt 

l2007HKuhlen & Faucher-Gieuerel 

120121 IRobertson et al.l 

20151). Once the ionizing photon 


production is determined we need to estimate its escape 
fraction fesc from galaxies to ionize the IGM. Direct ob¬ 
servational constraints of fesc are difficult, especially at 
z > 6 because of the high opacity of the intervening 
hydrogen residuals on the line of sight. Many studies 
and deep surveys were dedicated to the search of ion¬ 
izing continuum (Lyman continuum, LyG) escape from 
z < 4 galaxies. Very few detections were reported how¬ 
ever and show very low escape fractions relative to the 


UV radiation, of the order of few percent 

Shanlev et al. 

2006 

: llwata et al.l 12009 

: ISiana et alJl201C 

; INestor et al. 

2013 

). More recently. 

de Barros et al.l (201511 renorted 


a spectroscopic detection of Lyman continuum emission 
in a z ^ 3.2 galaxy with a relative escape fractiorFl of 
fesc 65%. 

In Figure 1111 we show our determination of the UV 
luminosity density together with the most recent re- 
sul ts in the literature as a fun ction of redshift i nclud ¬ 
ing IBouwens et alJ (l2015b[) andiFinkelstein et al.l (|2014[1 
at z = 4-8. IBouwens et al.l (|2015bll at z = 10, 

the relative escape fraction the ratio between the fraction of 
escaping Lyman continuum photons and the fraction of escaping 
photons at 1500 A (Itjteidei et ai.llS)01li 
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iMcLeod et alJ (|2015D and lOesch et alJ (I2014D at z ~ 9. 
While their values are based on the integration of the UV 
LF down to a magnitude limit between Mjjy = —18 and 
Muv = “17, our integration limit is two magnitudes 
fainter at Muy = —15. There is a significant d iffer¬ 
ence between the result s of iBouwens et al.l ()2ni5bl i and 
iFinkelstein et ahl (|2014li due to the fact that the UV LF 
in the form er study extends one mag nitude deeper than 
the latter. IFinkelstein et al.l (1201411 did not make use 
of the full IR data available in the HUDF, which also 
explains the larger uncertainties in their faint-end slope 
constraints. The UV luminosity density of galaxies at 
z ~ 7 determined in the present work is clearly larger 
than previous determinations, owing to a steep faint-end 
slope and a very faint integration limit. We can now as¬ 
sess whether this UV production is sufficient to ionize the 
IGM. We use the photon emission rate per unit cosmolog¬ 
ical comoving volume required to maintain reionizat ion 
at a given redshift determined in (iMadau et al.lll99^ : 

JV(*) = (1^)“ (6) 

where C 30 is the clumping factor Chu = (^hii)/(^hii)^ 
normalized to Chu = 30, and nnn is the mean comov¬ 
ing hydrogen density in the Universe. This is the mini¬ 
mum value for which the ionizing emission balances the 
recombination rate. We show in Fig. m this limit, con¬ 
verted to UV luminosity density (cf. IMadau et al.lil999l : 
IBolton fc Haehnef3l2007f) assuming an escape fraction of 
the ionizing radiation of fesc = 20% and three differ¬ 
ent values for the clumping factor. Assuming a standard 
value of CffTT = 3 for the clumping fac tor of the IGM 
(|Pawlik et al.ll200^ iFinlator et al.l[^12li our determina¬ 
tion of galaxy UV luminosity density at z ^ 7 is sufficient 
to maintain reionization. We also show that an escape 
fraction as low as fesc = 10 — 15% (green hatched re¬ 
gion) is already sufficient to ionize the IGM at z ~ 7. 
At the same redsh i ft, the ionizing emissivity constraints 
of IBouwens et al.l (|2015aj) are close to our UV density 
constraints. 

With current observations we are not able to put better 
constraints on the faint-end of the UV LF at z ~ 8, hence 
on the ionizing emissivity of galaxies at this redshift. The 
conclusions from the deep blank fields are still affected 
by large uncertainties that prevent any strong claims re¬ 
garding the con tribution of galaxies to the ionizing bud¬ 
get a t z ~ 8 (jBouwens et al.l l2015bl : IFinkelstein et al.l 
120141) . Future HST observations of the remaining HFF 
clusters might add better constrains on the UV luminos¬ 
ity density at z ~ 8 with additional highly magnified 
faint galaxies. 

6 . SUMMARY 

Combining HST observations of lensing clusters and 
parallel fields of the Hubble Frontier Fields program, we 
computed the galaxy UV luminosity function between 
z = 6 and z = 8. We assembled a large sample of about 
250 galaxy candidates in this redshift range using the Ly¬ 
man break photometric selection. In the cluster cores, we 
corrected the deep HST images for intra-cluster light and 
bright cluster galaxies light using a median filtering be¬ 
fore object detection, while performing the photometry 


in the original images. This technique helps the detec¬ 
tion of faint galaxies contaminated by cluster light but 
do not increase significantly the number of galaxies in the 
cluster center because of the very small volume probed 
in those regions that have very high magnification. 

• Using the latest lensing models produced by 
the CATS team fo r the three cl usters A2744 
(iJauzac et al.l I2015D . MACS0416 (iJauzac et al.l 
I2014D . and MACS0717 (Limousin et al. 2015), 
we have performed completeness simulations in 
the source plane that take into account all lens¬ 
ing effects. Thanks to the lensing magnification, 
the completeness function extends down to fainter 
magnitudes in the cluster fields than in the paral¬ 
lels. On the other hand, the total survey volume 
goes from ~ 4 arcmin^ in the parallel fields down 
to 0.6-1 arcmin^ in lensing fields. 

• We computed the UV luminosity function for all 
individual fields and for the combined sample at 
z ^ 7. The lensing magnification allow us to 
extend the LF down to an absolute UV magni¬ 
tude limit of Mjjy = —15.25, which is more than 
two magnitudes deeper than any study the deep 
blank fields. Most importantly, we show that 
the faint-end slope remains very steep at a = 
—2.041 q at such faint intrinsic luminosity. When 
excluding MACS0717 cluster, the uncertainties on 
the faint-end slope decrease to CTq ^ 0.1. The 
best Schechter fit yields a characteristic magnitude 
of M* = -20.891);;^^ and Log((;i*)=-3.54t°;^®. 
This is in good agreement with most of the re¬ 
cent results in blank fie lds (IBouwens et al.ll^OlSbl : 
IFinkelstein et ahl I2014D , which were limited to 
Muy = —17.5 and early results of lAtek et ^ 

(I 2 OI 5 I) . 

• Our determination of the UV LF at z ^ 8 does 
not reach beyond Muy = —18 because we detect 
only two galaxies in the two fainter magnitude bins. 
Albeit with large uncertainties due to small num¬ 
ber statistics and incompleteness uncertainties, the 
currently determined LF points confirm the results 
of previous studies. 

• On observational grounds, we integrate the UV LF 
down to a magnitude limit of Mjjy = —15 and 
find Log(p( 7 y) = 26.2 ± 0.13 erg s“^ Hz“^ Mpc“^. 
Assuming standard values for the ionizing conver¬ 
sion factor f^ion and the clumping factor Chu = 3, 
the ionizing budget of galaxies would be sufficient 
to maintain the IGM ionized by z ~ 7 provided 
the ionizing escape fraction from galaxies is greater 
than 10%. Future observations of the HFF pro¬ 
gram will allow us to decrease the uncertainties on 
the UV luminosity density at z ^ 7 and perhaps 
improve the constraints on the UV luminosity den¬ 
sity at z ^ 8, which are mostly based on blank field 
observations. 

With the help of gravitational lensing we have pro¬ 
duced the best constraints currently available on the UV 
luminosity function and the contribution of galaxies to 
the IGM reionization at z ~ 7. We clearly demonstrate 
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here the great potential and the feasibility of peering into 
the early universe through these cosmic telescopes. The 
remaining HFF clusters will certainly help improve even 
more our constraints on the UV LF at 2 : ~ 7 and in par¬ 
ticular at z ~ 8 to better assess the role of galaxies in the 
cosmic reionization process. Our results also show the 
great promise of future programs targeting lensing fields 
and pave the way to observing programs with the James 
Wehh Space Telescope or the Wide-Field Infrared Survey 
Telescope (WFIRST), which are scheduled for launch in 
the near future. 
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